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Abstract. The present article investigates the convergence of a 
class of space-time discretization schemes for the Cauchy prob- 
lem for linear parabolic stochastic partial differential equations 
(SPDEs) defined on the whole space. Sufficient conditions are 
given for accelerating the convergence of the scheme with respect 
to the spatial approximation to higher order accuracy by an appli- 
cation of Richardson's method. This work extends the results of 
Gyongy and Krylov [SI AM J. Math. Anal, 42 (2010), pp. 2275- 
2296] to schemes that discretize in time as well as space. 



1. Introduction 
For a fixed r G (0, 1), we consider the equation 

(1-1) = <-i + W«f + fi)r + J2 (Ml^vU + gl,) ef 

P =i 

for % G {1, . . . ,n} and (u),x) G Q x Gh with a given initial condition, 
where Gh is the space grid 

G h := {Ai/i + • • ■ + X p h; Ai, . . . , X p G A U (-A)} 

with mesh size h G R\{0} for a finite subset A C R d , for integer d > 1, 
containing the origin. For a fixed T G (0, oo) we define the time grid 

T T := {U = ir; i G {0, 1, . . . , n}, rn = T}, 

partitioning [0, T] with mesh size r, and note that v h = v h (u,t,x) 
depends on the parameter r as well as h, since we have used the 
convention of writing v\ in place of v h (ti) for G T T . In particu- 
lar, let = Au> p (tj_i) := w p (ti) — w p (tj_i) be the zth increment of 
w p with respect to T T , where, for integer d\ > 1, {w p ) d p =1 is a given 
sequence of independent Wiener processes carried by the stochastic 
basis J 7 , J-(t), P) that is complete with respect the filtration Fit) 
for t G [0, T\. For each i G {0, ...,n}, the L^ 1 and M i ' p are differ- 
ence operators given by L^<p := a^5h,\S-h^4> and Mf' p := b* p 5h,\<t>, 
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for p G {1, . . . ,d\}, where repeated indices indicate summation over 
X,/i G A. We assume that = a^(x) and b 4 A = (bf(a;))^i are real- 
valued V x ^-measurable functions on Q x T T x R d for all A, fi G A and 
further that = af A . Here P denotes the a- algebra of predictable 
subsets of Q x [0, oo) generated by Fit) and B = B(R, d ) denotes the 
a-algebra of Borel subsets of Tt d . The spatial differences above are 
defined by 

0(x + /iA) - 
d h>X (p{x) : = 

for A G R d \ {0} and by the identity for A = 0. We note that from this 
definition one can obtain both the so called "forward" and "backward" 
differences as h can be positive or negative. 
Together with (11. ip we consider 

di 

(1.2) Vi = + {dvi + f l )r + J2 {M^Vi-t + gU) 

P =i 

for i G {1, ...,n} and (u,x) G x R d with a given initial con- 
dition. Here Ci = C{tj) and Ai P = .M p (ti) are second order and 
first order differential operators given by C(t) := a Q/3 '(t)D a Dp and 
Ai p (t) := b ap (t)D a , respectively, where the summation is over at, /3 G 
{0, 1, . . . , d} and where D a = d/dx a , for a G {1, . . . , d}, while Do is 
the identity. For each a and we assume that a al3 (t) = a a "(t,x) and 
6 a (t) = (b ap (t, x))p =1 are real-valued V x ^-measurable functions on 
n x [0,T] x R d , and further that a a ^{t) = a^ a {t) for all t G [0,T]. 

Equations (II. II) and (II ,2p represent discrete schemes for approximat- 
ing the solution to the Cauchy problem for 

di 

(1.3) du(t,x) = (Cu(t,x) + f( y t,x))dt+^2(M p u( y t,x)+g p ( y t,x))dw p (t) 

for (u,t,x) G Q x [0,T] x R d with a given initial condition uo(x) = 
m(0, x). Under certain compatibility assumptions, equation (11.11) rep- 
resents an implicit space-time scheme for approximating the solution 
to the Cauchy problem for (II .3p by replacing the differential operators 
with finite differences and by carrying out an implicit Euler method in 
time. In a similar fashion, (jl.2p represents an implicit Euler method 
for approximating the solution to the Cauchy problem for (11.31) in time. 
Second order linear parabolic SPDE such as (11.31) arise in the nonlinear 
filtering of partially observable diffusion processes as the Zakai equa- 
tion ( [T3l [16j [221 [I])- Since analytic solutions to (II. 3p are difficult to 
obtain, there is a keen interest in providing accurate numerical schemes 
for its solution. 

Our aim is to show that the strong convergence of the spatial dis- 
cretization for the space-time scheme (II .ip to the solution of the Cauchy 
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problem for (II. 3p can be accelerated to any order of accuracy with re- 
spect to the computational effort. In general, the error of finite differ- 
ence approximations in the space variable for such equations is propor- 
tional to the mesh size h, for example, see [201 121] • We show the strong 
convergence of the solution of the space-time scheme to the solution of 
the time scheme (11.2)1 can be accelerated to higher order accuracy by 
taking suitable mixtures of approximations using different mesh sizes. 

This technique for obtaining higher order convergence, often referred 
to as Richardson's method after L.F. Richardson who used the idea to 
accelerate the convergence of finite difference schemes to determinis- 
tic partial differential equations (PDE) (see [T71 HB]), falls under a 
broadly applicable category of extrapolation techniques, for instance 
see the survey articles [21 E] ■ In particular, in [121 dH [TO] Richardson's 
method is implemented to accelerate the weak convergence of Euler 
approximations for stochastic differential equations. Recently, in [I] 
Gyongy and Krylov considered a semi-discrete scheme for solving (II. 3p 
which discretized via finite differences in the space variable, while al- 
lowing the scheme to vary continuously in time, and showed that the 
strong convergence of the spatial approximation can be accelerated by 
Richardson's method. The current paper extends these results to the 
implicit space-time scheme ( ll.ip . 

We must mention that for the present scheme one cannot also ac- 
celerate in time unless certain commutators of the differential operator 
A4 P in equation (II. 3p vanish, see [3|. For deterministic PDE we plan to 
address the simultaneous acceleration of the convergence of approxima- 
tions with respect to space and time in a future paper. Results concern- 
ing acceleration for monotone finite difference schemes for degenerate 
parabolic and elliptic PDE are given in [5], however our scheme is not 
necessarily monotone. 

In the next section, we present our assumptions as well as some 
preliminaries. Then in Section [3] we record the main results, namely 
Theorems 13.1) 13.2) and 13.3) the last of which says that the convergence 
of the spatial approximation can be accelerated to any order of accu- 
racy. In Section|4]we provide results which will be needed for the proofs 
of Theorems 13.11 and 13.21 In particular, we recall the solvability of the 
space-time scheme (II. ip . for the convenience of the reader, and present 
a new contribution — an estimate for the supremum of the solution to 
the scheme in appropriate spaces that is independent of h, the spatial 
mesh size. In Section Owe give the proof of a more general result and 
show that it implies Theorem 13.21 and hence Theorem 13.11 

We end with some notation that will be used throughout this work. 
Let £ 2 (Gh) be the set of real- valued functions (f> on G% such that 

MWH^E w*)i 2 <°° 

x€G h 
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and note that this notation will also be used for functions in i 2 {Rf). 

For a nonnegative integer m, let W™ = H^ n (R ) be the usual 
Hilbert-Sobolev space of functions on R d with norm || • || m . We note 
that for L 2 = L 2 (R d ) = the norm will be denoted by || • || - We use 
the notation D l (f> for the collection of all Zth order spatial derivatives 
of (p. Let 

W™(T) := L 2 (Sl x [0,T],V,W™) 

denote the space of W™- valued square integrable predictable processes 
on Q x [0, T]. These are the natural spaces in which to seek solutions 
to (HJ). 

2. Preliminaries and Assumptions 

We begin by setting some assumptions on our operators and recalling 
well known results concerning the solvability and rates of convergence 
for our schemes. In particular, we will discuss an t 2 {Gh) notion of 
solution and an L 2 notion of solution and recall an important lemma 
relating these function spaces. 

An L 2 -valued continuous process u = (u(t)) t< ={o t T] is called a general- 
ized solution to (jl.3p if u G W\ for almost every (u,t) G Q x [0,T], 




u{t)\\\dt < oo 



almost surely, and 

(u(t), 0) = [\(a°P - D a a a ?)D p u{s) + f(x), 0) - (a^Dpu, D a (j>) ds 
Jo 

+K 0) + V [\m p u(s) + go(s), 0) dw<>{s) 
P =i Jo 

holds for all t G [0,T] and <p e C%°(R d ). 

Assumption 2.1. For each {oj, t) G Q x [0, T] the functions a a/3 are m 
times and the functions b a are m + 1 times continuously differentiable 
in x. Moreover there exist constants K , . . . , K m+ i such that for I < m 

\D l a aP \ < Ki 

and for I < m + 1 

\D l b a \e 2 < Ki 

for all values of a, (3 G {0, . . . ,d} and (cu, t, x) G Q x [0, T] x H d . 
Assumption 2.2. There exists a positive constant k such that 

d 

(2a Q/3 - b ap b pp )z a z p > k\z\ 2 

a,/3=l 

for all (u,t,x) G Q x [0,T] x R d , z G K d , and p G {!,..., d x }. 
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Assumption 2.3. The initial condition uq G L 2 (£2, Tq, W™ +1 ), the 
space of W™ +1 -valued square integrable T^-measurable functions on Q. 
The f and g p , for p G {1, . . . , c?i}, are predictable processes onQx [0, T] 
taking values in W™ and W™ + , respectively. Moreover 

E [ (\\f(t)\\l+\\g(t)\\ 2 m+1 )dt + E\\u \\l +1 <oo, 
Jo 

where \\g{t)\\l-=T.U W^YWl 

Under Assumptions 12.11 12. 2\ and 12. 3[ the existence of a unique solu- 
tion u e w™ +2 (T) to is a classical result (see for example [151 H2] 
or Theorem 5.1 from [TT]). 

Remark 2.4. We note that by Sobolev's embedding of W™ C C&, the 
space of bounded continuous functions, for m > d/2 we can find a 
continuous function of x which is equal to uq almost everywhere for 
almost all u G Q. Likewise, for each (u,t) G Q x [0,T] there exists 
continuous functions of x which coincide with /(i) and g p {t) for almost 
every x G R d . Thus, if Assumption 12.31 holds with m > d/2 we assume 
that Uq, f{t), and g p {t) are continuous in x for all t G [0, T]. 

For a nonnegative integer m, let m := m V 1 and A := A \ {0}. We 
place the following additional requirements on our space-time scheme. 

Assumption 2.5. For all u G £1, for i G {0, . . . , n}, for X, p G A 0; 

and for v G A: the a Xfl are m times continuously differentiable in x; 
the a 0u and a u0 are m times continuously differentiable in x; and the 
b u are m times continuously differentiable in x. Moreover there exist 
constants A , . . . , such that for A, p G A and j < m we have 

\D j a x> *\ < Aj 

and for A G A and j < m we have 

\D j a X0 \ < Aj, \D j a 0X \ < A j} and \D j b x \ < Aj 

for all (oj, x) G x H d for i G {0, . . . , n}. 

Assumption 2.6. There exists a positive constant k such that 

(2a v - b xp b pp )z x z^ >kJ24 
A,/^eAo AeAo 

for all (u, x) G Q x R^, i G {0, . . . , n}, p G {1, . . . , d\}, and numbers 
z x , A G A . 

For (11. ip to be consistent with (11.31) we also require the following. 

Assumption 2.7. For i G {0, . . . , n} 

„oo ^oo 

E„A0\a | „°Ai, „a0 , „0ct 

AeA AteA 
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£ ofAV = of, 

A,/i6Ao 

and 

aga 

for all a, j3 G {1, . . . , d} and p G {1, . . . , d{\. 

Remark 2.8. If Ao is a basis for and Assumption 12.71 holds then 
Assumption 12.11 implies 12.51 and 12.21 implies 12.61 with m = to. 

A solution v h = (t>f )™ =1 to (11.11) with an ^(G/J-valued J-o- m easurable 
initial condition Vq is understood as a sequence of £ 2 (Gh)-valued ran- 
dom variables satisfying (11.11) on the grid Gh- The following result is 
well known and we provide it for the sake of completeness. 

Theorem 2.9. Let f and g p be T%- adapted £ 2 (Gh) -valued processes and 
let Vq be an T§-measurable l 2 {Gh) -valued initial condition. If Assump- 
tion \MJk holds then (11.11) admits a unique £ 2 (Gh) -valued solution for 
sufficiently small t. 

Proof. By Assumption 12.51 for each % e {1, . . . , n}, equation (II. ip is a 
recursion with bounded linear operators on £ 2 (Gh). In particular, for 
each h the operator norm of tL h is smaller than a constant less than 
1 for sufficiently small r, independently of w 6 O. Hence (I — tL k ) 
is invertible in £ 2 (Gh) for sufficiently small r, by the invertibility of 
operators in a neighborhood of the (invertible) identity operator /. 
Therefore, for i e {l,...,n} we are guaranteed an £ 2 (G/ l )-valued 
satisfying {I — rL^)<f> = if) for all if> G £ 2 (Gh) and moreover this solution 
is easily seen to be unique. Thus we can construct a unique solution 
to the scheme iteratively. □ 

The rate of convergence of the solution v h of (11.11) (and v of (II. 2p ) 
to the solution u of (11.31) with initial condition u is known. In [6], [71 
8], Gyongy and Millet obtained the rate of convergence for a class of 
equations in the nonlinear setting of which our schemes are a special 
case. Namely, in the situation of Remark 12. 8| if Assumptions 12.11 12.21 
and 12.31 hold with a a @, b a , f, and g p all Holder continuous in time with 
exponent 1/2 then 

£max V V \6 h>x (v^(x)- Ui {x))\ 2 h d 

|A|<m+l x£G h 
n 

+ EtJ2 E ^\8hAvfa)-Ui(x))\ 2 h d <N(h 2 +T) 

1=1 \\\<m+2x£G h 

for sufficiently small r, h G (0, 1), and for a constant that is indepen- 
dent of h and r. The principal interest of this paper is to investigate 
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higher order convergence with respect to the spatial discretization, that 
is, to obtain an estimate, similar to the above, with a higher power of 
h by applying Richardson's method. 

While it is natural to seek solutions to (II. ip on the grid, carrying 
out our analysis on the whole space will have certain advantages when 
it comes to providing estimates for solutions to our schemes. Indeed, 
we observe that (11. ip is well defined not only on Gh but for all x G R d . 
Therefore, we introduce an alternate notion of solution. A solution to 
fll.ip on Q x T T x R d with an L 2 -valued J-"o-measurable initial condition 
Vq is a sequence v h = (t>^)™ =1 of L 2 -valued random variables satisfying 
fll.ip . In a similar spirit, solutions to f ll.2p with the appropriate initial 
condition are understood as sequences of H^-valued random variables 
satisfying fll.2p in W^ 1 . The next result follows immediately from the 
considerations in the proof of Theorem 12.91 

Theorem 2.10. Let f and g p be Ti-adapted L 2 -valued processes and 
let Vq be an J^Q-measurable I? -valued initial condition. If Assumption 
\2.5\ holds then fll.ip admits a unique L 2 -valued solution for sufficiently 
small t. 

By Sobolev's embedding theorem, for / > d/2 there exists a linear 
operator I : W\ — > Cb such that <p(x) = Icj){x) for almost every x G H d 
and sup xgR d \I(j)(x)\ < N\\<p\\i for all G W% where N is a constant. 
We recall the following useful embedding of W\ C £ 2 (Gh) from [I]. 

Lemma 2.11. Let I > d/2 and \h\ G (0,1). For all (f) G W\ the 

embedding 

(2.i) \mx)\ 2 \h\ d < mm 

x€G h 

holds for a constant N that depends only on d and I. 

Proof. For z G R d let B r fx) := {x G R d ; \x — z\ < r}. By the embed- 
ding of W\ into Cb, for G Cb we have 

\4>(z)\ 2 < sup <j) 2 (z + hx) 
xe-Bi(o) 

<iVV/i 2|a| / \{D a cf>)(z + hx)\ 2 dx 



\a\<l 

< \h\ 2lahd [ \(D a (j))(x)\ 2 dx 

\a\<l ^ B h( z ) 

<N\h\~ d Y) I \(D a <l>)(x)\ 2 dx 

\ n \<l jB h{») 



(0) 



\a\<l 

for a constant N depending only on d and I and thus 

M(p h) = E \<l>(z)\ 2 \h\ d < N J2 E / \(D a Mx)\ 2 dx 

JB h (z) 



z£G h \a\<lz€G h 
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which yields the desired embedding. □ 

We will show that the restriction of a continuous modification of 
an L 2 -valued solution to (11. ip to the grid Gh is also a solution in the 
£ 2 (Gh) sense. Thus we will carry out our analysis in the whole space 
and obtain estimates independent of h in appropriate Sobolev spaces 
for the L 2 -valued solutions of (11.11) and (11.21) . 

We provide the aforementioned Sobolev space estimates in Section 
IH We then use these estimates in Section \5\ to prove the main results, 
which are the focus of the next section. 



3. Main Results 

To accelerate the convergence of the spatial approximation by Richard- 
son's method we must have an expansion for the solution v h to (II. ip 
with initial data Vq = uq in powers of the mesh size h. This relies on 
the possibility of proving the existence of sequences of random fields 
v(°\x), v^(x), v( k \x), for x G R d and integer k > 0, satisfying 
certain properties. Namely, v^°\ . . . are independent of h; is 
the solution of (11.21) with initial value uq] and an expansion 

(3.1) *) = tT^W + ^ 

i=o J - 

holds almost surely for i G {1, . . . ,n} and x G Gh-, where R T,h is an 
C-2 (Gh) -valued adapted process such that 

(3.2) Emax sup \R T - h {x)\ 2 < Nh 2(k+1) lC m 

x£G h 



i<n 



for 



n 



K m := E||u ||m+1 + Er Y)(\\fiWm + ll^llm+l) < 00 



lm+1 ^ LJI / j \ 
i=0 

and a constant independent of r and h. 

Our first result concerns the existence of such an expansion. 

Theorem 3.1. // Assumptions [2J\, El ELS HOE and\2l\ hold 
with 

d 

m = m>k + l + - 

for an integer k > then expansion (13. ip and estimate ( 13. 2 p hold for 
a constant N depending only on d, d\, A, m, K , . . . , K m+ \, A , . . . , 
A m , k, and T . 

In the proof of Theorem 13.11 as v is defined not only on Gh but for 
all x G R d , we will see that one can replace Gh in (13.21) with R d . We 
also note that in the situation of Remark 12.81 if Assumptions 12.11 and 
12.21 hold with m > k + 1 + d/2 then the conditions of Theorem 13. II are 
satisfied. 
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Taking differences of expansion (13. ip clearly yields 

Sh^(x) = ]T ^5 hA v?\x) + 5 h , x Rl' h (x) 

for any A = (Ai, ...,A P ) G A p , for integer p > 0, where A := {0} 
and 6h,x '■= $h,Xi x ■ • • x S^x . The bound on 5h,\R T,h is not obvious, 
nevertheless we have the following generalization of the above theorem. 



Theorem 3.2. If the assumptions of Theorem \3.1\ hold with 

d 

m = m>k + p + l + - 
for a nonnegative integer p then for A G A p expansion (13.11) and 
£max sup \6 h x R? h (x)\ 2 + Emax\h\ d V \5 htX R? h (x)\ 2 < Nh 2 ^K m 

hold for a constant N depending only on d, d\, K, m, K , . . . , K m+ \, 
A , . . . , A m , k, and T. 

The proof of Theorems 13. II and 13.21 appear in Section [5] following the 
considerations in the next section. Currently we shall discuss how to 
implement Richardson's method to obtain higher order convergence in 
the spatial approximation, extending the result from [I] to the space- 
time scheme. 

Fix an integer k > and let 



(3-3) 

j=0 



2~ih 



where v 2 solves, with 2 in place of h, the space-time scheme 
(II. ip with initial condition Uq. Here (3 is given by (f3 , /3i, . . . , (3k) '■— 
(1, 0, ... , O)^ 1 where V' 1 is the inverse of the Vandermonde matrix 
with entries V ij := 2~^- l ^ j -^ for i, j G {1, . . . , k + 1}. Recall that v (0) 
is the solution to (II. 2p with initial condition u . 

Theorem 3.3. Under the assumptions of Theorem \3.1[ 

(3.4) £max sup \v?(x) - vf\x)\ 2 < N\hf k+1) K, m 

l < n xeG h 

for a constant N depending only on d, d 1; A, m, K , . . . , K m+1 , A , 

.... . Kj . (XTid T . 



Proof. By Theorem 13.11 we have the expansion 

v 2-'h = v (o) + y ^0) + r ,2-*h h k+i 



1=1 
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for each j G {0, 1, . . . , k} where r T ' 2 ~ Jh := h~^ k+ ^R T ' 2 ^ h . Then 



v 



o=o / j=0 i=l ' j=0 

i=l ' j=0 j=0 

= u (0) + S£ /3jr ,2-ih h k+l 
3=0 

since X^=o ft = 1 anc ^ Sj=o ft2 _lJ = for each i G {1,2,..., k} by the 
definition of (/?o, • • • , At). Now using the bound on i? T ' h from Theorem 
13. II together with this last calculation yields the desired result. □ 

One can also construct rapidly converging approximations of deriva- 
tives of v(°\ That is, if the conditions of Theorem 13.11 hold instead 
with 

d 

m = m>k + p + l + - 

for nonnegative integers k and p then Theorem 13.31 holds with 5h 1 \V h 
and 8h,\v^ in place of v h and v^°\ respectively, for A G A p . Therefore, 
using suitable linear combinations of finite differences of v h one can 
construct rapidly converging approximations for the derivatives of v^°\ 
In the next section, we present material that will be used to prove 
the main results in this section. In particular, we provide estimates 
for the L 2 -valued solutions of (11.11) and fll.2p in appropriate Sobolev 
spaces. 

4. Auxiliary Results 

We include the following bound, which is given for the continuous 
time case in [4], for the convenience of the reader. 

Lemma 4.1. If Assumptions [Ql and \2.b\ hold then for all <fi G L 2 

Qi{<l>):= / 2<P(x)L^(x) + J2\M^ P m\ 2 dx 
J *> d P =i 

AeA 

for all i G {1, . . . ,n} and for a constant N depending only on k, A , 
A\, and the cardinality of A. 

Proof. First observe that for /i G A the conjugate operator in L2 to 
8-h,ix is —8h,n- Notice also that 
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where Th :fl ip(x) = ift(x + hfi). Thus by simple calculations Q = + 

Q(2) + Q(3) + Q(4) where 



Q?\<P) :=- f £ ^ ~ $ P K P )M)KMx) dx, 
QS 3) (0) := 2 / (o°V (s) + <j>{x) £ K A %,a0 + a° A <U»(*)) dx, 



and 



Qf } (0) := / (b?V(:r) + 2 £ b^<M M 0(*)) 



>h.\<»\\o 



By Assumption 12. 6j 

Qi 1 ^) <-«£>* 

AeA 

and by Assumption 12.51 Young's inequality, and the shift invariance of 
Lebesgue measure, 

Q?\*)<^\MM + nm\\1 

AeAo 

for each j G {2,3,4} with a constant A depending only on the cardi- 
nality of A, k, A and, for j = 2, also on A\. □ 

We also recall the following discrete Gronwall lemma. Note that we 
use the convention that summation over an empty set is zero. 

Lemma 4.2. For constants K G (0, 1) and C , if (aj)™ =0 is a nonnega- 
tive sequence such that aj < C+K J2i=i a « holds for each j G {0, . . . , n} 
then aj < C(l - K)~ j for je{0,...,n}. 

Proof Let bj — C + KJ2l=i k and note that (1 — K)bj = bj-\. Then 
aj < bj for j < n by induction, since a < C = b and 

i-i j-i 
^•(l-A) <C + ^2a i <C + ^2b i = b j {l-K), 

i=l i=l 

assuming a.,_i < Therefore a, < bj = C(l — A)~ J for each j < n 
and for AT G (0,1). ' □ 

The following provides a Sobolev space estimate for solutions to the 
space-time scheme that is independent of h. For an integer m > 0, 
denote by W™(r) the space of W™- valued predictable processes on 
Q x T T such that 



i=i 
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and note that we write 



n d± 
V \^ II „P||2 



\g\ m = Er^^\\g 
i=i p=\ 



for functions g = (g p )pL 1 . 

Theorem 4.3. For fx E A and p E {1, . . . , dx}, let ^ e W^(r) . 
// Assumption holds then for each nonzero h there exists a unique 
solution v E W™(r) of 

(4-1) = Vi-x + j>s + f?> + E( M "-i^-i + ^i)c p 

^teA p=i 

/or an?/ W^^-va/wed To-measurable initial condition uq. Further, if 
Assumption \2.6i is also satisfied then 

n 

£max \\viWn + Er /] y~] II^M^IIm < ^^ll^ollm+i 

i<n ' * ' ' 

(4.2) i=1 XeA 



+NErj2(\\m 2 m + n 



i=0 



holds for a constant N that depends only on d, dx, m, A, Aq, . . . , A„, 
k, and T . 

Proof. By Theorem 12.101 the existence of a unique sequence of L 2 - 
valued random variables solving (11. ip is known. For f tl ,g p E W™(r) 
and an W™ + - valued initial condition z/ , there exists a unique sequence 
of W™- valued random variables satisfying (14. ip . The estimate (14. 2 p can 
be achieved easily with a constant iV depending on h, so in particular 
the solution is in W™(r). 

Next we prove the estimate (14. 2 p for a constant independent of h. 
For convenience we denote .ft™ := r SILodl/illm + llftllm)- Considering 
equalities of the from a 2 + b 2 = 2a(a — b) — \a — b\ 2 we note that (14. ip 
implies 

IHIo ~ = 2(Vi,Vi - Vi-x) - \\Vi - Vi-x\\l 

= 2(i/<, lS + ft)r + 2(zvi, M&u^ + gUtf 
+ 2(ui - Ui-i, M^Vi-i + g P i-x)t,i u i-i\\o 

= 2{v h L\vi + ft)r + 2(1*.!, M&Vi-i + gl^ 

+ \\(M^x + gl^Wl - + f?\\y 

where here and in what follows we suppress the sums over p E A and 
p E {1, . . . , dx}. Summing up over i, we have 



(4.3) \\ujWl < HWl + Uj+Tj + Jj 



J 
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where 



2(1/,, Lfi/, + /nT 
X, :=^2(^ 1 ,M^> l _ 1 + ^_ 1 )e 



and 



~pi|2 

3 j Ho- 



^:=£ll(^i-i + *£-i)tf 
i=l 

By an application of Ito's formula, it is easy to see that for rr,p £ 

^+1^+1 = (Aw"(U))(Aw»(U)) = YZ - Yr + t5 7TP 
for all i G {1, . . . , n) where 



TT"{t):= / (w*(s)-w; {s) )dwO(s)+ / (w'(s)-w p is) )dw*(s), 
Jo Jo 

j(s) is the piecewise defined function taking value j(s) = i for s G 
[it, (i + 1)t), and 8 wp = 1 when it = p and otherwise. Thus can write 

Jj = Jf ] + Jf ] where 

:=EEn^-i+Ciii^ 

i=i p=i 

and 

7T,P=1 

Then we note that since Lemma 147X1 holds for all t G [0, T], in particular 



«i + < ^E \\ s vto\\l + ^E (g<("0 + k ft) + ^hi 



_ , H 

AeA i=i 



< Cr||^b||? + Nt E H|g - |r E E H^llo + ^ 

i=l i=l AeA 

where here || folio < Cll^oll? an d iV is a positive constant depending 
only on k, A , Ax, and the cardinality of A. Thus we can replace 
equation (14. 3p by 

j j 

mi +^EE n^iio < nt\Wo\\i + ntJ2 Nio 

(.4-4) «=i AeA i=i 

+iV^+X i + J r J (2) 
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for a constant N that depends only on k, A , Ai, C, and the cardinality 
of A. 

Next we observe that 



j 

EX i = tl 2 J nd E { E (^-i(M^-i + gU)% | Fi-x)) dx = 



since is independent of T{ and Vi, M i ,p i/ i , and g p are all Ti- 

measurable for % G {0, . . . ,n}. Similarly, we see that Ejj 2) = since 
the expectation of the stochastic integral is zero. Therefore, taking the 
expectation of f)4.4j) we have that 



E\\ Vj \\l + ErJ2J2 II V^llo < NEt\\v \\1 + NE% 

j 

+NET^2pi\\ 



(4.5) 



8=1 



for each j 6 {1, . . . ,n}. Excluding for the time being the difference 
term on the left hand side of (14.51) and applying Lemma 14.21 we obtain 



EWujWl < N(ET\\v \\l + ES§)(l-NT)- j 
and, since (1 - Nr)~ j = (1 - N^)~ j < (1 - N-)~ n < C'e NT , we have 



(4.6) max£7||i/<||jj < NEr\\u \\l + iVE^ 



for a constant iV here that depends only on the parameters k, A , Ai, 
T, and the cardinality of A. Using equation (14. 6[) we can eliminate the 
last term on the right hand side of (14. 5p . In particular, we have 



E ^J2J2W 5h ' xUi Wo <NEt\\v \\1 + NES%. 

i=l AeA 
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Using the Davis inequality we can bound max and max|X|. 

Namely, 

Emax\j!> 2) \ 

i<n J 

di , r T ^ 1/2 

<3J2 E \ II^t(.) + ^)II2IK^7W + ^.)IIS^ p >W 

w,p=l ^ Jo 

di , pT ^ 1/2 

7r,p=l WO ' 



T,P= 

di 



<CYE I max || AfVi/i + # 

7T,p=l ^ 



X 



- / K(W) + <( S )II>^ s )-<mI 



where C is a constant independent of r and h that is allowed to change 
from one instance to the next. Therefore, 



di 

P||2 

i Ho 



£max|j; (2) | < diCj" rE m&xWM^Vi + g; 

i<n ' J i<n 

(4.7) 

+7 E ^ ril<w^) + < w lloK(^)-<( s) l 2 ^ 

by Young's inequality. The first term on the right hand side of ( 14. 7p is 
bounded from above by the sum over alH e {1, . . . , n}, hence 



£ r£ max ||JW* "V, + II 2 , < Er £ ||M^ + 

p=l ' i=0 



< CEr^J^ IK,a*4o t- ll.'Min'- 

i=0 AeA 

and the second term on the right hand side of (14. 7p yields 

d\ r T 



7T,p=l J ° 

7T,P=1 \ \ J0 



s) lloK( s ) ~ W l(s)\ dslfys) 



^ Er E(EH^llo + ll^llo) 

i=Q AeA 
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by the tower property for conditional expectations. Combining these 
estimates we see that Emax\J'^\ is estimated by terms already ap- 
pearing on the right hand side of (I4.4p for a constant N that depends 
also on d\. 

Similarly, we note that 

3 

X j = ]T 2(1*-!, M&Vi-! + ^ti)A<! 
i=l 

ti 

T h <P „ , , _U nP \ JnnPt 



2 / (i/ 7(a) , M 7( ' s p )Z / 7(s) + ^ (s) ) dw'(s) 

«/ 



Applying the Davis inequality 



EmaxlXil < gJ^eI / ll^lloll^M^rM + 0?( s )llo ds 
- p =i < J o 

di I ( f T 1 1 ^ 

- 6 E E (^< a „ x II Ui II { J o II M i& u -r(') + 1 1 o rfs j 

and then Young's inequality 

1 " 
(4.8) £max|2i| < -^max + C^r V(V ||5 M ^||g + 

i=Q AeA 

we see that max |X| is also estimated by terms already appearing on 
the right and side of (14.41) . 

Returning to (I4.4p and taking the maximum followed by the expec- 
tation we have 

n 

EmaxpiWl + ErS^y^ \\6 h xUi\\l < NEt\\v \\\ + NER™, 

i<n ' J ' J 

i=l AeA 

using (14. 6p and the estimates on E max | ,7^1 and £7max|X|. Thus 
holds when m = 0. 
If m > 1 we differentiate (14.11) with respect to x l and introduce the 
notation <fi for the derivative of a function <p in the direction x l for 
I £ {1, . . . , d}. Then (14.11) becomes 



-i + ^ S h,\S- htll i>i + ft)r 
A,peA 

(4-9) ^ 

p=l AeA 

where / M := / M for nonzero /i, f° := f° + a XfJ, Sh,\5^h,fj,v and g p := 
g p + b Ap 5/ lj Ai / - Recalling that <9„ = for G A , we proceed as 
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before but now using the inequality 



i=l i=l 

n N n 

< eErJ2\\ D Kx^\\l + J Er Yl 

i=l i=l 



r. II 2 



* 110 



which holds for arbitrary e > and N depending only on |/i|. This 
leads to the following 



ElfaW* + EtJ2J2 ll^llo < NEr\\u \\l + NE&™ 

(4.10) 1=1 AeA 

3 3 

+ 2d ET IP^IIo + NEtJ^MI 

i=l AeA i=l 

for each ar, Z G {1, . . . , d}. Summing up over each direction x , the term 
with factor l/2d can be seen to be estimated by other terms already 
appearing on the right hand side of f)4.10p . Then by the same procedure 
as before we obtain 

n 

(4.11) Emax\\Dui\\l + ET V V ||L>^ A ^||o < NEt\\u \\1 + NE&l 

i<n ' J ' J 

i=l AeA 

which proves the theorem when m = 1. 

Assuming that m > 2 and that (14.21) holds for each integer p < 
m in place of m, then we can differentiate (14. ip (p + 1) times and, 
repurposing the notation for the (p+ l)th order derivatives of </> with 
respect to x, we obtain (14.91) with different f° and g p . Namely, the 
f° will be the sum of f° and linear combinations of certain zth order 
derivatives of a Xfl together with certain (p + 1 — i)th order derivatives 
of 5h,\5-h !f iV, for integer i < (p + 1). As before, the L 2 -norms of the 
(p + 1 — z)th derivatives of S^xS-h,^ are dominated by the L 2 -norms 
of the (p + 2 — i)th derivatives of 5h,\ v which are in turn less than the 
W2 +1 -norm of S^x 1 *- After similar changes are made in g p we obtain 
the counterpart of (14. lip which then yields (14. 2 p with (p + 1) in place 
ofm. □ 



We also have the following Sobolev space estimate for solutions to 
the implicit time scheme. Let m be a nonnegative integer. 

Theorem 4.4. Let f e W^(r) and g p G W™ +1 (r). If Assumptions 
\2.1\ and \2.2\ hold then (II. 2p has a unique solution v G W™ +2 (r) for a 
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given W™ +1 -valued J-'o -measurable initial condition Vq. Moreover 

n 

^maxlHlm+l + Et Y\ INIm+2 < NE \\vofm+i 
l<1l ' 

1=1 

n 

+NETj2(\\fi\\ 2 m + \\9i\\ 2 m+ i) 

i=0 

holds for a constant N depending only on d, di, m, Kq, . . . , K m+ i, k, 
and T . 

Proof. Proving the solvability of (11 .2p reduces to solving the elliptic 
problem 

di 

(/ - r&)vi = Vi-x +rfi + J2 ZKMUvi-t + gU) 

for each i G {1, . . . , n} where I is the identity. That is, we claim that 
A := (/ — tC) is a W™- valued operator on W™ +2 such that Ai is 

(i) bounded, i.e. ||^4j0||^ < fC||0||^ t+2 for a constant K, 

(ii) and coercive, i.e. (Ai4>,(p) > A||0||^ +2 for a constant A > 0, 

for every i G {1, . . . , n} and for all <f> G W™ +2 , where (•, •) denotes the 
duality pairing between W™ +2 and W™ based on the inner product 
in W™ + . Then by the separability of W™ +2 , there exists a countable 
dense subset {ej}JL 1 such that for fixed p > 1, <f> p = Y7j=i c j e j f° r 
constants Cj where 4> p G W™ +2 is not identically zero. We fix i and for 
every ip G W™ consider A acting on (f) p , that is (A(f> p , ef) = ef) for 
all j G {1, . . . ,p}. Taking linear combinations we obtain {A(f> p , 4> p ) = 
(ip, cf) p ) from which we derive 

by the coercivity of A and an application of the Cauchy-Schwarz in- 
equality. Thus ||0 p || m +2 < xll^llm and hence (by the reflexivity of 
W™ +2 ), there exists a subsequence pk such that <\> Ph converges weakly 
to (p and in particular (Ai4> Pk , ef) — > (A<p, ef) for every j. Therefore for 
every ip G W™ there exists a (f> G W / 2 m+2 satisfying Ai4> = ip for every 
i G {1, . . . , n}. Moreover, this solution is easily seen to be unique. 

Using the existence and uniqueness to the elliptic problem in each 
interval, we note that 

di 

by Assumption 12.31 and therefore there exists a v i G W™ +2 satisfying 

di 

(/ - T&fa =V + Tfx + ^2( M V + 9 P )ei- 

p=0 
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Further, assuming that there exists a V{ e W™ +2 satisfying (II. 2p we 
have that 

di 

Vi + rf i+1 + Y^Mivt + g^Ui G 

by the induction hypothesis and Assumption I2.3| and therefore there 
exists a v i+1 e W™ +2 satisfying (11.21) . Hence we obtain t> = (fj)" =1 
such that each v { G H^ 2 m+2 satisfies (|T72|) . 

It only remains to prove the claim concerning ellipticity of A. By 
Assumption 12.11 clearly A% is a bounded linear operator for each i. We 
see that 

{A<j>A)^{i4>A)-r{C4>A) 

= U\\l +1 -r(C<PA) 
where, by Assumptions 12.11 and 12. 2| 

(£0, 0) = ((a ^ - D a a^)D p <l>, 0) - (a^D^, D a <f>) 

< cu\\ 2 m+1 - ^U\\ 2 m+2 

in the W™ +1 inner product for a, [3 <E {I, . . . , d} and for a constant C 
depending on K$ and K\. Therefore 

(M, 0) > %rM\\l+2 + (1 - ^)il0ll^ + i > \ r||0||^ +2 

for sufficiently small r and hence (ii) is satisfied. 

To prove the estimate, we use a method similar to that in the proof 
of Theorem 14.31 to arrive at (14. 3p with v in place of u, £ in place of L h , 
Ai p in place of M h ' p , all in the W^" 1- -norm instead of the L 2 -norm. 
Again, we decompose J into and using the processes Y np (t) 
that arise by applying the Ito formula to the product of increments of 
the independent Wiener processes. However this time, instead of using 
Lemma |4JJ we observe that 



j j 3 

^+^ (1) <^r^||^|| 2 rt+1 -|r^||^|| 2 „ +2 +iVr^(||/ J || 2 n +||^ 



2 ' / j II ~* nm-t-z . - ■ - / J \i|.' ' inn \\a% II m+1 / 
i=l i=l i=0 



since 



di 

2v{x)Cv{x) + \M p v(x)\ 2 dx <(e- kC)\\v\\1 +2 + C\\v\\ 2 m+1 

Rd T^i 

for e > by the considerations above. Therefore we have that 

j 

2 i ^V^||.,.||2 / A7IL, 112 , Apr. , Ar^(2) 



\ V j\\m+l 



+ rJ2 INIm+2 < N\\vo\\ 2 m+1 + NT, + NJ 3 



i=l 

J 



2 

«llm+l/ 

i=0 
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and the estimate follows by considering the maximum and then tak- 
ing the expectation. Moreover, with the estimate, it is clear that the 
solution v e W™ +2 (r). □ 

We will use the theorem above to obtain estimates in appropriate 
Sobolev spaces for a system of time discretized equations. For i e 
{0, . . . , n} and an integer p > 1, let 

AeA 

and let 

A,/iSAo 3=0 AeA 

+(p+ir i E«^ +1 > 

AteAo 

-MP := (P+I)- 1 ^ 

aga 

iMp) ._ rU) 



and 



i=o J 



<^(p)p . = M h,p _ '_LM U)p 
1 1 7 1 

i=0 J 

where A p j is defined by 

(-i)p-3 

(4.12) A M - 



(j + l)\(p-j + l)V 

For p > 1, the values of C^cj) and A4^ p (j) are obtain by formally taking 
the pth derivatives in h of L h <j) and M h,p cf) at /t = 0. 

For a positive integer k < m, the sequences of random fields t> W, . . . , 
nee( i ec [ j n f|3.i|) will be the embeddings of random variables taking 
values in certain Sobolev spaces obtained as solutions to a system of 
time discretized SPDE. Namely, as the solutions to 

(4-13) 



i=i 

for p G {1, . . . , k} where C = p(p — 1) • • • (p — I + is the binomial 
coefficient and z/°) is the solution to (11. 2p from Theorem H] 
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Theorem 4.5. Let Assumptions \2.1\ \2.2\ \2.3[ and \2.5\ hold with m = 
m > k > 1 and let z/ 0) G W™ +2 (r) be the solution to (ll.2p with initial 
condition uq from Theorem \4-4\ Then the system (I4.13P with initial 
condition 



m+2— p. 



T . 



has a unique set of solutions (z/^)p =1 such that each v&> G W 
Moreover for each p G {1, . . . , /c}, 

^max||^||^ p + Er^||^||^ +2 _ p 
(4-14) „ 1=1 

<NErJ2(\\m 2 m + H\\ 2 m+ i) 

i=0 

holds for a constant N depending only on d, d\, A, m, Kq, . . . , K m+ i, 
A , . . . , A m , k, and T . 



Proof. For convenience let 



v 



and 



p 



3=1 

where we write = G {p)p . 

Observe that for each p G {1, . . . , k) the equation for z» in (jUSD 
depends only on v"> for I < p and does not involve any of the unknown 
processes with indices I > p. Therefore we shall prove the solvabil- 
ity of the system and the desired properties on i/W recursively using 
Theorem I4.4L 

For p = 1, we have 

(4.15) v p = !/« + (Az/W + if y + + Gft)^ 

P =i 

Since G W™ +2 (r), for m = m we have that F« G W^fV) 
and G^ 1 ' G W™(r) by Assumption 12.51 Hence by Theorem 14.41 there 
exists a unique v^ 1 ' G W™ +1 (r) satisfying (14.151) with initial condition 
z/q 1 '' = 0. Further, is estimated by (I4.14p and thus Theorem 14.51 
holds with p — 1. 

Now we assume that for m > k > 2 and p G {2, . . . , k} we have 
unique . . . , z/( p_1 ) solving (14.131) for = . . . u^" 1 ' = with the 
desired properties. In particular, observe that for j G {1, . . . ,p} 

(4.16) [£^z>^)] m _ p < N ] m+2 - (p-j) 
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and for each p G {1, . . . , d{\ 

(4.17) [M^u^U-p+.i < 

for a constant N. Therefore it follows that G W™~ p (r) and G 
W™ _p+1 (r). Applying Theorem 14.41 yields the existence of a unique 
solution z/( p ) G W™~ p+2 (r) that satisfies (14.131) with initial condition 
Uq = 0. Together with (I4.16P and (14.171) the estimate from Theorem 
14.41 implies that (14. 14j) holds. Further, the uniqueness of each 
follows from Theorem 14.41 □ 

For the convenience of the reader we record the following lemma and 
two remarks from [3J that will be used in proving the error estimates. 

Lemma 4.6. Let (ft G W 2 +1 and if) G W 2 +2 for a nonnegative integer 
p and let X, /x G A . Set 



Then we have 

(4.18) 

and 

QP 



QP 



and d 



e p d p x +l c 



{x + hex) d9 



(4.19) 



(any 



-1 pi 



[ [ (e 1 d x -e 2 d fl ) p d Xt ,ip(x + h(e 1 x-e 2 p))de 1 de 2 

Jo Jo 



for almost all x G R d for each h G R. Furthermore, for integer I > 
if (ft E W p+2+l and if) G W p+W th 



en 



(4.20) 

and 
(4.21) 



h? 



< 



p+1 



(p + 2)! 



0~h,\0~-h,fi'f 



< N\h\ p+1 



l+p+3i 



i=0 j=0 

where A+j is defined by (14.121) and N depends on X, fi, d, and p. 

Proof. It suffices to prove the lemma for (ft, if) G C^°(Ti, d ). For p = 0, 
formula f 14 . 1 8 j) is obtained by applying the Newton-Leibniz formula to 
(ft(x + 6hX) clS db function of 9 G [0, 1]. Namely, 



(x + hX) — <ft(x) 



u=x+h\ 



Dj(ft(u) du 



h 



X j D j( ft(x + OhX) dd 



e=o 



and therefore Sh,\(ft(x) = d\(ft(x + OhX) d6. Applying the Newton- 
Leibniz formula again yields (14. 19[) with p = 0. After that, for p > 1 



ACCELERATED SPATIAL APPROXIMATIONS 23 

one obtains (I4.18P and (I4.19P by differentiating both parts of these 
equations written with p = 1. 

Next by Taylor's formula for smooth f(h) we have 

v-^ W d? 1 & d p+1 

f(h) = Y———f(0) + - (h-9) p ———f(9)d9. 

Applying this to 

S h ,x(f>(x) = [ d X (f>(x + 9hX) dd 
Jo 

as a function of h we see that 

w 

77 

3=0 



^ x ) = E(jTT)!^ +l0(x) 

p// (l-9 2 ) p 9 P+1 d p x +2 4>(x + h9 l 9 2 X)d9 1 d9 2 . 
• Jo Jo 



h P+i rWi 

+ 



Now to prove ()4.20j) . it remains only to use that by Minkowski's integral 
inequality the W^-norm °f the last term is less than the W^-noim. of 



d p (f) times 



\h\p+i /■! /-i i \h\ p+1 



pi 



Similarly, by observing that the value at h = of the right hand side 
of (1419]) is 

p«E^ +1 fir i+ V(x), 

j=0 

we see that the left hand side of (14.211) is the W^-norm of 
i.p+1 ri ri ri 

—T (l-03) p (Oidx-9 2 d ll ) p+1 d x ^(x+h9 3 (9 1 \-9 2 fi))d9 1 d9 2 d9 3 , 
P- Jo Jo Jo 

which yields ffl~2Tj) . □ 

For integers / > and r > 1, denote by Wjf 2 the Hilbert space of 
functions on R d such that 

(4.22) U\\l r>h := £ \\5 hM x .-. x ^011? < oo 

Ai,...,A r eA 

and set W l ff 2 = W\. Then for any <p 6 H^ +r we nave 



where N depends only on |A | 2 := 2~^AeA l^| 2 an( ^ r - 
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Remark 4.7. Formula (I4.18P with p = and Minkowski's integral 
inequality imply that 

||^,a0||o < ||<9a0||o- 

By applying this inequality to finite differences of <fi and using induction 
we can conclude that W l 2 +r C W l { 2 . 

Remark 4.8. Owing to Assumption 12. 7\ for i G {0, . . . ,n} we have 
that Cf^ = Ci and M.f^ p = M.^. Also by Lemma FPjI and Assumptions 
[23] and EH for G W% +2+l and if; G W% +W we have 

\\O h ^\\i < N\hr l U\\ l+p+ s 

and 

ll^^ll, < N\h\*> +1 \\<f>\\ l+p+2 

for a constant N depending only on p, d, I, A , . . . , A h and A. 

For integers k, I > 0, let z/°), i/W, . . . , i/W be the functions from 
Theorem 14.51 We define 

fc 1,7 

(4.23) rj> h := if - - £ 



J'=l J 

for i G {1, . . . ,n} where v h is the unique L 2 - valued solution to (11. ip 
that exists by Theorem 14.31 with initial condition uq, data f° — f and 
/" = 0,/i6A O . 

Lemma 4.9. Let Assumptions \2.1\ \2.2\ \2.3[ and \2.5\ hold with tn = 
m = I + k + 1 for integers k, I > and /et r r,/l be defined as in equation 
flQ3|) . Taen rj' fc = ; r T > h G W™~ fc (r) and 

/or i £ {l,...,n} where 
and 



i=o J 



and, moreover, F T ' h G W|(r) and G rAp G W^fr). 

Proof. By Theorem 14.31 the solution to the space-time scheme z/ h G 
W™(t) and by Theorem 14.41 the solution to the time scheme G 
W™ +2 (r). Therefore r T ' h G W?(t) when jfe = and, by Theorem 1431 
r r,h e w^- fe (r) when jfc > 1. 
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Observe that 

' //' //•' x ^ //' x ^ //' ^ „, 



Z — / j| Z.^ 7! Z — / ^| Z — < j\ 

i=0 ' j=0 J ' i=0 ' j=l J ' 

=yy Kt^^ 

Z.^ Z.^ o\ 7 1 

i=l j=0 ' J ' 

^ ^ - 1)\ 

i=l 3=1 yj > 



7 U - 7 
i=l j=l J v J ' 

where summations over empty sets are zero. Therefore, we can rewrite 

F T ' h as 



Similarly, observe that 



7 ^ 7 



y " y " M a)P U (i) = yy " =: j^p 

i=0 j=0 J i=l 3=1 v J/ 

and therefore 



A- k ■ 



G T 'V = M h ' p u w - .M V 0) + V —M h ' p v ij) - y —M p u {j) - J rAp . 

Z — / 7 1 Z — <• 7 1 

Thus, following from Remark 14.81 and Theorem 14. 5[ F T,h G W^r) and 
Gr,fc,P e Wf, +1 (r). " □ 

With the previous considerations, we are now prepared to prove the 
main results. 

5. Proof of Main Results 

We prove a slightly more general result which implies Theorem 13.21 
Here we suppose that m = m. 

Theorem 5.1. Let Assumptions {Ml I£M {KM {KM {KM and^7J\ hold 
with m = I + k + 1 for integers l,k > 0. Then for r T ' h as defined in 
(05j) we have 

n 

(5.1) Sn^llr^llf + Sr^^ll^'l? < N\h\ 2 ^}C m , 

i<n ^— ' ^— ' 

i=l AeA 

where N depends only on d, d\, A, m, K , . . . , K m+ i, A , . . . , A m , k, 
and T . 
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Proof. Recall that by Lemma 14.91 we have that F T,h G W^fY) and 
qtAp e w 1 2 +1 (t). Then the left-hand-side of ([EE]) is dominated by 



"AP||2\ 

1/ J 



(5.2) NEr^m^Wi + W^ 

i=l 

due to Lemma 1431 and Theorem 14.31 To estimate (I5.2p we observe that 
for j < k, by Remark 14.81 we have that 

\\lSJy\n>\ |Pj \\l +k -j+3 — iV|Al| ||Z/j \\m+2-j, 

and combining this result with Theorem 14.51 yields 

n 

ErJ2\\ F P h \\i < N\h\ 2{ - k ^K, m . 
i=i 

The bound on G T,h ' p can be obtained in a similar fashion, yielding the 
desired result. □ 

Now set R T,h := Ir T,h where / is the embedding operator from 
Lemma 12.111 We have the following corollary to Theorem 15.11 which 
implies Theorem 13.21 



Corollary 5.2. If the assumptions of Theorem ] 5.1\ hold with I > p+d/2 
for a nonnegative integer p then for A G A p 

Emax sup \8 hjX Rl' h {x)\ 2 < Nh 2{k+l) lC m 



and 

Emax V \6 htX R-' h (x)\ 2 \h\ d < Nh 2 ^ k+1) K, m 

i<n ^— ' ' 

xeG h 

hold for a constant N depending only on d, d\, A, m, Kq, . . . , K m+ x, 
A , . . . , A m , k, and T. 

Proof. Using Sobolev's embedding of W^ p into Cf, and Remark 14.71 
Theorem 15.11 implies 

Emax sup \8h\Rl' (x)\ 2 < CE max \\r1 ,h \\ 2 _ h 

< C'Emax\\rT> h \\? 

i<n 

< Nh ^)jc m 

where C and C are constants depending only on m and d, and N is 
a constant depending only on m, c?i, K, A, iT , . . . , iT m+1 , and T. 
Similarly, by Lemma [2.111 above and Remark 14.71 

£max V \8 h>x Rl' h {x)\ 2 \h\ d < C£max ||<5 M i^||f _ 

i<n ' ' i<n 

x&G h 



< C'E max \\rj' h \\ 2 



i<n 



< Nh 2 ^ k+1) ]C r 
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□ 

For the / : W\ — > Cb from Lemma 12.111 Theorem 13.21 follows by 
considering the embeddings v h := Iu h , where v h is the unique L 2 - 
valued solution to (II. ip with initial condition Uq, and 

v (j) . = jj/y) 

for 

j G {0, . . . , fc}, where z/°) is the unique L 2 - valued solution to (II. 2p with 
initial condition u and the processes v^ l \ . . . , are the solutions to 
the system of time discretized SPDE (I4.13P as given in Theorem 14. 51 By 
Theorem l4.3| u h is .^-adapted and Jy2 -vame d for all z G {1, . . . , n}. For 
each j G {1, . . . , k} the i/^) are iy p+1+fc - valued processes by Theorem 
14.51 Since I > d/2 and p+ 1 — k> d/2 the processes v h and i;^ are well 
defined and clearly (14.231) implies (13. ip with v h in place of v . That 
is, we have the expansion for a continuous version of the L 2 -valued 
solution. 

To see that Theorem 13.21 indeed follows from Corollary 15.21 we must 
show that the restriction of the L 2 - valued solution to the grid Gh, a 
set of Lebesgue measure zero, is indeed equal almost surely to the 
unique £ 2 (G/ l )-valued solution that one would naturally obtain from 
(II. ip . That is, we must show that 

(5.3) $(x) = «?(*) 

almost surely for all i G {1, . . . , n} and for each x G Gh where v h is the 
unique ^-adapted £ 2 (G ? ? i )-valued solution of (ll.ip from Theorem 12.91 
Therefore, for a compactly supported nonnegative smooth function <fi 
on R d with unit integral and for a fixed x G Gh we define 

My) ■= (^) 

for y G R d and e > 0. Recall, by Remark 12. 4| that we can obtain 
versions of uq, f, and g p that are continuous in x. Since v h is a L 2 - 
valued solution of (II. ip for each e, almost surely 

/ v?(y)Mv)dv= I Vi-i(y)4>e{y)dy + r I (LW + fi)(y)<l> e (y)dy 

jR d jR d JR d 

for each i G {1, . . . , n}. Letting e — > 0, we see that both sides converge 
for all z G {1, . . . , n} and u> G Q. Therefore almost surely 

= «?-!(*) + (LW(x) + / i (x))r + Y<( M t p ivlAx) + 

p=l 

for all i G {l,...,n}. Moreover by Lemma 12.1 1[ the restriction of 
v h , the continuous version of v h , onto Gh is an £ 2 (G/ l )-valued pro- 
cess. Hence ( 15. 3 p holds, due to the uniqueness of the £ 2 (Gh)-valued 
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J-i-adapted solution of (II. ip for any £ 2 (G/ l )-valued J-o- m easurable ini- 
tial data. This finishes the proof of Theorem 13.21 

We end with the following generalization of Theorem 13.31 

Theorem 5.3. If the assumptions of Theorem \3.B hold with p = and 
v h as defined in (\3.3\i then 

Em&x sup \v^(x) — vf\x)\ 2 

+ £max V \vHx)-vf\x)\ 2 \h\ d < N\h\ 2{k+1) IC m 

i<n ^— ' 

x&G h 

for a constant N depending only on d, d\, A, m, K , . . . , K m+1 , A\, 
. . . , A m , k, and T . 

This follows from Theorem 15.11 and the definition of v h . 
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